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Abstract 

An effect generated by the nonexponential behavior of the survival 
amplitude of an unstable state in the long time region is considered. 
In 1957 Khalfin proved that this amplitude tends to zero as t — > oo 
more slowly than any exponential function of t. For a time-dependent 
decay rate j(t) Khalfin's result means that this j(t) is not a constant 
for large t but that it tends to zero as t — > oo. We find that a similar 
conclusion can be drawn for the instantaneous energy of the unsta- 
ble state for a large class of models of unstable states: This energy 
tends to the minimal energy of the system £ m in as t — > oo which is 
much smaller than the energy of this state for t of the order of the 
lifetime of the considered state. Analyzing the transition time region 
between exponential and non-exponential form of the survival ampli- 
tude we find that the instantaneous energy of a considered unstable 
state can take large values, much larger than the energy of this state 
for t from the exponential time region. Taking into account results 
obtained for a model considered, it is hypothesized that this purely 
quantum mechanical effect may be responsible for the properties of 
broad resonances such as a meson as well as having astrophysical and 
cosmological consequences. 
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1 Introduction 



Within the quantum theory the state vector at time t, | $(£)), for the physical 
system under consideration which initially (at t — to = 0) was in the state 
|$) can be found by solving the Schodinger equation 

iK^\$(t)) = H\$(t)), |$(0)) = |$), (1) 

where |3>) € TC, 7i is the Hilbert space of states of the considered sys- 

tem, || |$(i))|| = || |$)|| = 1 and H denotes the total selfadjoint Hamiltonian 
for the system. If one considers an unstable state |<&) = \<p) of the system 
then using the solution \<p(t)) of Eq. for the initial condition |0(O)) = \<p) 
one can determine the decay law, V</>(t) of this state decaying in vacuum 

V^t) = \a{t)\\ (2) 

where a(t) is probability amplitude of finding the system at the time t in the 
initial state |</>) prepared at time to = 0, 

a(t) = tfMt)). (3) 

We have 

a(0) = I. (4) 

From basic principles of quantum theory it is known that the amplitude 
a(t), and thus the decay law V^it) of the unstable state |0), are completely 
determined by the density of the energy distribution u{£) for the system in 
this state [T] 

a{t) = [ u{S) e~i Et d£. (5) 

JSpec.{H) 

where u>(£) > 0. 

Note that ((SJ) and (j4]) mean that there must be 

o(0) = / uj(£)d£ = 1. (6) 

From the last property and from the Riemann-Lebesgue Lemma it follows 
that the amplitude a(t), being the Fourier transform of u(£) (see (jSJ)), must 
tend to zero as t — > oo PQ. 
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In [2] assuming that the spectrum of H must be bounded from below, 
(Spec.(H) > — oo), and using the Paley- Wiener Theorem [3] it was proved 
that in the case of unstable states there must be 

\a(t)\ > Ae~ ht \ (7) 

for \t\ — ► oo. Here A > 0, b > and < q < 1. This means that the decay 
law V^>{t) of unstable states decaying in the vacuum, (jSJ), can not be described 
by an exponential function of time t if time t is suitably long, t — > oo, and 
that for these lengths of time V^(t) tends to zero as t —>■ oo more slowly 
than any exponential function of t. The analysis of the models of the decay 

processes shows that V<t>{t) — exp( — £-), (where 7° is the decay rate of the 
state |0)), to a very high accuracy for a wide time range t: From t suitably 
later than some T ~ £0 = but T > to up to t ^> = \ and smaller than 

t — t as , where t as denotes the time t for which the nonexponential deviations 
of a(t) begin to dominate (see eg., [2], [1] [8]. From this analysis it follows 
that in the general case the decay law V^(t) takes the inverse power-like 
form t~ x , (where A > 0), for suitably large t > t as ^> r^. This effect is in 
agreement with the general result (j7j). Effects of this type are sometimes 
called the "Khalfin effect" [9]. 

The problem of how to detect possible deviations from the exponential 
form of Vcj,{t) at the long time region has been attracting the attention of 
physicists since the first theoretical predictions of such an effect. [El QUI ED] 
Many tests of the decay law performed some time ago did not indicate any 
deviations from the exponential form of V<f,(t) at the long time region. Nev- 
ertheless, conditions leading to the nonexponetial behavior of the amplitude 
a(t) at long times were studied theoretically [12] - [16]. Conclusions follow- 
ing from these studies were applied successfully in an experiment described 
in the Rothe paper [T7], where the experimental evidence of deviations from 
the the exponential decay law at long times was reported. This result gives 
rise to another problem which now becomes important: If (and how) devia- 
tions from the exponential decay law at long times affect the energy of the 
unstable state and its decay rate at this time region. 

Note that in fact the amplitude a(t) contains information about the decay 
law V<f,(t) of the state \<f>), that is about the decay rate 7° of this state, as 
well as the energy £® of the system in this state. This information can be 
extracted from a(t). Indeed, if \<p) is an unstable (a quasi-stationary) state 
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then 

a(t) = e~tt S * ~ (8) 

So, there is 

b * 2^~ in dt a(ty [J) 

in the case of quasi-stationary states. 

The standard interpretation and understanding of the quantum theory 
and the related construction of our measuring devices are such that detect- 
ing the energy £® and decay rate 7° one is sure that the amplitude a(t) has 
the form (jSj) and thus that the relation (Q occurs. Taking the above into ac- 
count one can define the "effective Hamiltonian" , h^,, for the one-dimensional 
subspace of states 7i\\ spanned by the normalized vector |0) as follows (see, 
eg- P]), 

, def .f. da(t) 1 , . 

h * = lh -dT^- (10) 

In general, can depend on time t, = h<f,(t). One meets this effective 
Hamiltonian when one starts with the Schrodinger Equation ([T]) for the total 
state space TC and looks for the rigorous evolution equation for the distin- 
guished subspace of states TC\\ C TC. In the case of one-dimensional H\\ this 
rigorous Schrodinger-like evolution equation has the following form for the 
initial condition a(0) = 1 (see, [18] and references one finds therein), 

fcMfi = hS)a{t). (11) 

Relations (TlOT) and (flTJ) establish a direct connection between the amplitude 
a(t) for the state |</>) and the exact effective Hamiltonian h<f,(t) governing the 
time evolution in the one-dimensional subspace Tl» 3 \<p). Thus, the use of 
the evolution equation pip or the relation fflOl) is one of the most effective 
tools for the accurate analysis of the early- as well as the long-time properties 
of the energy and decay rate of a given quasi-stationary state \4>{t)). 

So let us assume that we know the amplitude a(t). Then starting with 
this a(t) and using the expression (TlOl) one can calculate the effective Hamil- 
tonian hcj)(t) in a general case for every t. Thus, one finds the following 
expressions for the energy and the decay rate of the system in the state |</>) 
under considerations, to be more precise for the instantaneous energy and 
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the instantaneous decay rate, (for details see: [T9]). 



8* = £ 6 (t) = ^(h 6 (t), (12) 
= 7^) = -29f(M*)> ( 13 ) 

where dt (z) and Q (z) denote the real and imaginary parts of z respectively. 

The deviations of the decay law V^it) from the exponential form can be 
described equivalently using time-dependent decay rate (IT3J) . In terms of 
such 70 {t) the Khalfin observation that V<p{t) must tend to zero as t — > oo 
more slowly than any exponential function means that J,j>(t) <C 7° for £ ^> t as 
and lim^oo ^(t) = 0. 

Using (flOl) and (fl^]) . (fT31) one can find that 



^(0) = (14) 

£*(*~t*) ^ 5; ^ ^(0), (15) 

7^(0) = 0, (16) 

7*(*~r*) ~ 7 J. (17) 

The aim of this talk is to discuss the long time behaviour of £<j,{t) using 
a(t) calculated for the given density u(8). We show that £<f,(t) — > 8 min > —00 



as t — > 00 for the model considered and that a wide class of models has similar 
long time properties: 8 ( j } {t)\ t ^ oo 7^ £?. It seems that, in contrast to the stan- 
dard Khalfin effect [2], in the case of the quasi-stationary states belonging 
to the same class as excited atomic levels, these long time properties of the 
instantaneous energy £$(€) have a chance to be detected, eg., by analyzing 
the properties of the high energy cosmic rays or the spectra of very distant 
stars. 

2 The model 

Let us assume that 

Spec.(H) = [8 m i n , 00), (where, £ m i n > —00), and let us choose u>(8) as 
follows 

N Is 

u(£) = uj bw (8) = — 9{£ - £ min ) — — , (18) 

2tt ( g _ £ ay + (?i.y 
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where iV is a normalization constant and 

6{£) - 



1 for £ > 0, 
for £ < 0, 



For such ojbw{£) using (jSJ) one has 



a(t) = — -5— e » d£, 

2vr A_ (£_£0 )2 + ( ^ )2 



where 



1 - 1 r ^ 



^ 2 W frain (£_£0 )2 + ( | )2 



Formula (IT9l) leads to the result 

-V 

i 



— x 

2tt 



a(t) = Ne h( £ i 1 2*)* |l 
-v9t 

x e~ £?i ^ - - £ TO i„ + - 

-^(-^S-^n-^)*)]}, 

where E\(x) denotes the integral-exponential function [20| [21]. 
In general one has 

a(t) = a exp (t) + a non (t), 

where 

a exp (t) = Ne~^ £ l ~ a non (t) = a(t) - a exp {t). 

Making use of the asymptotic expansion of E\ (x) [21] 



where | argz| < §7r, one finds 



N —iP.ff. .. 7$ ft 



- g wiw ) T ° /fix a + 

I £m.ir). I \~t / 



where fcj = £j - §7$, and 



2 '0' 

da(t) 1 



^ a(t) 



t— >oo 

_ .ft _ £° ~ £mm /ft\ 2 

£"min I \t J 

for the considered case (TTHT) of ujbw(£) (f° r details see [19J). From 
follows that 

K(V*)UJ = *~(*) 

_ ^0 — ^min /ft\ 2 

~ r — / 

— '-'vn.i.rt. ^ 



\h° — £ • I 2 \t 

| "(A '-"mm | x 17 



"mini 



where £™(t) = ^(*)|t-oo, and 

ft 



t— »oo 



The property (1261) means that 

*(M*)LJ = W < 

For different states \<p) = \j), (j — 1,2, 3, . . .) one has 

whereas in general 7° 7^ 7°. 

Note that from fl24l) one obtains 



/V 2 1 7 ° + t5 

+— sin[(£$-£ min )*] e "3T* _^ - 

, iV 2 (7^) 2 ft 2 



47T 2 J ft ? £min | ^ t 2 



Relations (1211) — (1251) become important for times t > t as , where t as denotes 
the time t at which contributions to l^)^^ 2 from the first exponential 
component in (l30l) and from the third component proportional to 4- are 
comparable, that is (see (1221) ). 

|a exp (t)| 2 ~ |a non (t)| 2 (31) 

for t — >■ oo. So t as can be be found by considering the following relation 

e -| t„ J_ (7^) 2 ^_ (32) 

Assuming that the right hand side is equal to the left hand side in the above 
relation one gets a transcendental equation. Exact solutions of such an equa- 
tion can be expressed by means of the Lambert W function [22] . An asymp- 
totic solution of the equation obtained from the relation (1321) is relatively easy 
to find [23]. The very approximate asymptotic solution, t as , of this equation 
for (-$) > 10 (in general for (-$■) — > oo) has the form 



J A, Las . w <--/> 



8,28 + 4 In 



+ 2 1n[8,28 + 4 1n( ^ ~ /"*" )] + ... . (33) 

3 Some generalizations 

To complete the analysis performed in the previous Section let us consider a 
more general case of u>(£) and a(t). For a start, let us consider a relatively 

simple case when lim£_ ) .£ min+ = f u; > and u)(£)\ £<£ = 0. Let 

derivatives uj( k '(£), (k = 0, 1, 2, ... , n), be continuous in [£ m i n , 00), (that is 
let for £ > £ m i n all u>^ k '(£) be continuous and all the limits lini£^£ m . n+ u^ k '{£) 
exist) and let all these u^ k '(£) be absolutely integrable functions then (see 

a(t) ~ -^e-^g^!)*^)*^ (34) 



k=0 



where 4 = f lim^£ mi „+ u/ fc) (£). 
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Let us now consider a more complicated form of the density uj{£). Namely 
let uj{£) be of the form 

lo {£) = {£- £ min ) x rj(£) e Li(-oo,oo), (35) 

where < A < 1 and it is assumed that r/ fc )(£), (A; = 0,1,..., n), 
exist and they are continuous in [£ m i n , oo), and limits lim £ ^ £min+ r]^(£) ex- 
ist, lim^oo {£ — £ rn i n ) x rj( k \£) = for all above mentioned k and 
= 0, then 



ih , —if . f 

\ . fc ^ III I II 1 



a(t) ~ — — A e a n (t) + ( - — ) a n ^(t) 



t 

t— >oo 



/j 



+ (- -T- J a„_ 2 (t) 

+ (- j) 3 a n _ 3 (t) + • • • > (30) 



where (compare [25l [26] ) 

n— fe — 1 



™r([+A) _,^±^) (!+*)/ V* ^ 

<*«-*(*) = 2^ /j e 2 % (t) , (37) 

r(z) is the Gamma Function and ?7q"'' ) = lim£_,.£ min+ r]^(£), r]^(£) = rj(£) 
and j = 0, 1, . . . , n. 

The asymptotic form of h u (t) for t —>■ oo for the a(i) given by the relation 
looks as follows 

^tt CO = ^ttCO It— = ^min — i "7 

-^(y) 2 + .... (38) 

In the more general case of a(i) (see, e.g. (1361) ) after some algebra the 
asymptotic approximation of a(t) can be written as follows 



■t N 

a(t) ~ e^^^T-^, (39) 



t^oo k=0 
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where £ > and are complex numbers. 
From the relation (l3T?j) one concludes that 

da(t) —itf ( % „ v-^ ,x cu 



- h E min {_«^_£ (e + ife) _g_}. (40) 



at 

t^oo k=0 



Now let us take into account the relation ( TTTT) . From this relation and 
relations p?]) , pOI) it follows that 

. . _ g?i do d"i , , , 

M*) ~ £ mm + J + J + J + ■ ■ ■ , (41) 

where di, d<i, e^, . . . are complex numbers with negative or positive real and 
imaginary parts. This means that in the case of the asymptotic approxima- 
tion to a(t) of the form the following property holds, 

lim hf{t) = £ mm < £j. (42) 

It seems to be important that results (jUJ) and (|42p coincide with the 
results ( 1251) — ( 1291) obtained for the density cJeiy(£) given by the formula 
( {TBI) . This means that general conclusion obtained for the other oj{£) defining 
unstable states should be similar to those following from (|25|) — (|29|) . 



4 Final remarks. 

Long time properties of the survival probability |a(t)| 2 and instantaneous 
energy £<f>{t) are relatively easy to find analytically for times t ^> t as even 
in the general it was shown in previous Section and [24J. It is much 

more difficult to analyze these properties analytically in the transition time 
region where t ~ t as - It can be done numerically for given models. 

The model considered in Sec. 2 and defined by the density u>bw(£), (|T8|) . 
allows one to find numerically the decay curves and the instantaneous energy 
e^it) as a function of time t. The results presented in this Section have been 
obtained assuming for simplicity that the minimal energy £ m i n appearing in 
the formula ffTS]) is equal to zero, £ m in — 0. So, all numerical calculations 
were performed for the density lobw{£) given by the following formula 

£W(£) = ^ 0(£) 7 * -o , (43) 
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Figure 1: Survival probability V^(t) = \a(t)\ 2 in the transition time region. 
The case % = 10. 

£° 

for some chosen Performing calculations particular attention was paid to 

the form of the probability |a(t)| 2 , i. e. of the decay curve, and of the instan- 
taneous energy e^{t) for times t belonging to the most interesting transition 
time-region between exponential and nonexponential parts of |a(£)| 2 , where 
the following relation corresponding with fl3Tl) and fl32|) takes place, 

\a exp {t)\ 2 ~ |a no „(t)| 2 , (44) 

where a exp (t), a non {i) are defined by fl22|) . Results are presented graphically 
below. 

A similar form of a decay curves one meets for a very large class of models 
defined by energy densities u>(£) of the following type (see [Sj 127]). 

u{£) = —0(£ - £ min ) (£ - £ min f '* f(£), (45) 

= ^Bw{£) {£ — £min) X f{£)i 

where A > 0, f(£) is a form-factor — it is a smooth function going to 
zero as £ — > oo and it has no threshold and no pole. The asymptotical 
large time behavior of a(t) is due to the term {£ — £ m in) x and the choice 
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17 18 19 20 21 l~<f> 



Figure 4: Survival probability V^{t) = \a(t)\ 2 in the transition time region. 



£2 



The 




Figure 5: Survival probability V${t) = \a(t)\ 2 in the transition time region. 
The case % = 100. 
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29.5 30.0 30.5 31.0 31.5 32.0 T <t> 



Figure 6: Survival probability V^t) = \a(t)\ 2 in the transition time region. 
The case 5 = 100. 




Figure 7: Instantaneous energy in the transition time region. The case 



'* - 100. 
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of A (see Sec. 3). The density u(£) defined by the relation (T4"5l) fulfills all 
physical requirements and it leads to the decay curves having a very similar 
form at transition times region to the decay curves presented above . The 
characteristic feature of all these decay curves is the presence of sharp and 
frequent oscillations at the transition times region (see Figs (JTJ, (j2J), (J4j) , 
©, © )• This means that derivatives of the amplitude a(t) may reach 
extremely large values for some times from the transition time region and 
the modulus of these derivatives is much larger than the modulus of a(t), 
which is very small for these times. This explains why in this time region the 
real and imaginary parts of h^(t) = £<j,(t) — |7^,(£), which can be expressed 
by the relation (1101) . ie. by a large derivative of a(t) divided by a very small 
a(t), reach values much larger than the energy £® of the the unstable state 

measured at times for which the decay curve has the exponential form. For 

e° 

the model considered we found that, eg. for = 10 and 5t^ < t < QOr^ 

the maximal value of the instantaneous energy equals £</,(£) = 89, 2209 £° 
and £ ( f > (t) reaches this value for t = t mXy 10 = 53, 94 and then the survival 
probability V^it) is of order V^,(t mXtW ) ~ 10~ 9 . 

The question is whether and where this effect can manifest itself. There 
are two possibilities to observe the above long time properties of unstable 
states: The first one is that one should analyze properties of unstable states 
having not too long values of t as . The second one is finding a possibility to 
observe a suitably large number of events, i.e. unstable particles, created by 
the same source. 

The problem with understanding the properties of broad resonances in 
the scalar sector (a meson problem [28]) discussed in [221130], where the hy- 
pothesis was formulated that this problem could be connected with properties 
of the decay amplitude in the transition time region, seems to be possible 
manifestations of this effect and this problem refers to the first possibility 
mentioned above. There is the problem with determining the mass of broad 
resonances. The measured range of possible mass of a meson is very wide, 
400 - 1200 MeV. So one can not exclude the possibility that the masses of 
some u mesons are measured for times of order their lifetime and some of 
them for times where their instantaneous energy £ a {t) is much larger. This 
is exactly the case presented in Fig. ([3]) and Fig. ([7]). For broad mesons the 
ratio -J is relatively small and thus the time t as when the above discussed 
effect occurs appears to be not too long. 

Astrophysical and cosmological processes in which extremely huge num- 
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bers of unstable particles are created seem to be another possibility for the 
above discussed effect to become manifest. The probability V${t) = |a(t)| 2 
that an unstable particle, say 0, survives up to time t ~ t as is extremely 
small. Let V<j>{t) be 

nWL ias ~ 10^, (46) 
where k ^> 1, then there is a chance to observe some of particles <p survived 
at t ~ t as only if there is a source creating these particles in N<f, number such 
that 

m)\t~ taa ^ » i. (47) 

So if a source exists that creates a flux containing 

~ 10 1 , (48) 

unstable particles and / ^> A; then the probability theory states that the 
number N surv unstable particles 

N surv = V^t)\^ tas M <t) ~ lO 1 "* > 0, (49) 

has to survive up to time t ~ t as . Sources creating such numbers of unstable 
particles are known from cosmology and astrophysics. The Big Bang is the 
obvious example of such a source. Some other examples include processes 
taking place in galactic nuclei (galactic cores) and inside stars, etc. 

So let us assume that we have an astrophysical source creating a suffi- 
ciently large number of unstable particles in unit of time and emitting a flux 
of these particles and that this flux is constant or slowly varying in time. 
Consider as an example a flux of neutrons. From (1331) it follows that for the 
neutron t™ s ~ (250r„, — 300r„), where r n ~ 886 [s]. If the energies of these 
neutrons are of order 30 x 10 17 [eV] then during time t ~ t™ s they can reach 
a distance d n ~ 25000 [ly], that is the distance of about a half of the Milky 
Way radius. Now if in a unit ot time a suitably large number of neutrons M n 
of the energies mentioned is created by this source then in the distance d n 
from the source a number of spherically symmetric space areas (halos) sur- 
rounding the source, where neutron instantaneous energies £ n {t) are much 

2 

larger than their rest energy £® = , mn Z „ , (m? is the neutron rest mass 

and v n denotes its velocity) have to appear (see Fig. (jSJ)). Of course this 
conclusion holds also for other unstable particles <p a produced by this source. 
Every kind of particles 4> a has its own halos located at distances df. a , 

4 a -^Jt\ (k = i,2,...), 
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Figure 8: Halos surrounding a source of unstable particles 
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Figure 9: Halos generated by a few not too distant sources S a , Sp. 



from the source. Radiuses d^ a of these halos are determined by the particles' 
velocities v^ a and by times tf' s k when instantaneous energies £$ a (t) have local 
maxima. 

Unstable particles Q forming these halos and having instantaneous en- 
ergies £tf, a (t) 3> £® a = m^c 2 have to interact gravitationally with objects 
outside of these halos as particles of masses m^ a (t) = j% £<f> a (t) ^> mj a . 
The possible observable effects depend on the astrophysical source of these 
particles considered. 

If the halos are formed by unstable particles emitted as a result of internal 
star processes then in the case of very young stars cosmic dust and gases 
should be attracted by these halos as a result of a gravity attraction. So, 
the halos should be a places where the dust and gases condensate. Thus in 
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the case of very young stars one may consider the halos as the places where 
planets are born. On the other hand in the case of much older stars a presence 
of halos should manifest itself in tiny changes of velocities and accelerations 
of object moving in the considered planetary star system relating to those 
calculated without taking into account of the halos presence. 

If the halos are formed by unstable particles emitted by a galaxy core 
and these particles are such that the ratio -fir is suitably large inside the 

halos, then rotational velocities of stars rounding the galaxy center outside 
the halos should differ from those calculated without taking into account the 
halos. Thus the halos may affect the form of rotation curves of galaxies. (Of 
course, we do do not assume that the sole factor affecting the form of the 
rotation curves are these halos). Another possible effect is that the velocities 
of particles crossing these galactic halos should slightly vary in time due to 
gravitational interactions, i. e. they should gain some acceleration. This 
should cause charged particles to emit electromagnetic radiation when they 
cross the halo. 

Note that the above mentioned effects seems to be possible to examine. 
All these effect are the simple consequence of the fact that the instantaneous 
energy £^ a (t) of unstable particles becomes large compared with £® a and for 
some times even extremely large. On the other hand this property of £^ a (t) 
results from the rigorous analysis of properties of the quantum mechanical 
survival probability a(t) (see fl3]) ) and from the assumption that the energy 
spectrum is bounded from below. 
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